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Abstract 



We consider the factorization method in the spontaneously broken gauge theories such 
as the electroweak theory. 

Using the off-shell lU-boson density matrix formalism we demonstrate that the factor- 
ization conditions are completely under control. 

The main point of this paper is the presence of the "interference" or "crossed" terms in 
lUlU-scattering process which exhibits itself in the dependence on the relative azimuthal 
angle between two scattering planes. 

To illustrate our general consideration we use well-known example — a single Higgs 
boson production. The origine of the quantitative failure of the WW -effective method, as 
originally propossed, is given and it is shown how to use it correctly. 
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I. Introduction 



In order to search for a new physics at future e + e~-colliders or at LHC we must well 
understand the fusion mechanism especially via vector W, Z bosons. 

Our paper is devoted to the applicability of the effective H^-boson approximation (EWA)- 
method. This approximate method was proposed for the calculation of WW(ZZ) -fusion 
processes, for which there were no hopes to find exact analytical answers in general case. 

Originally [1, 2] only longitudinally polarized H^-boson beams were taken into account 
by EWA. It was found that the result of calculation of the cross-section for Higgs-boson 
production via H^H^-fusion by EWA (in the region of W-boson energy E w > m H > my/) 
is much larger than the exact total cross-section [3]. The contribution of the transversaly 
polarized W-heams to this process only increases the discrepancy between this classical term 
and the exact total cross-section. 

A lot of efforts were spent to resolve this discrepancy between exact computer calculations 
and the EWA-method. But the results were not very succesful - all the corrections were 
found to be small and the EWA-method was announced to be correct only for heavy Higgs- 
boson (M H 3> m w ) and only in high energy limit E w 3> M H [4, 5]. 

In this paper we suggest a resolution of this problem based on the advanced version of 
EWA-method developed in our paper [6]. Actually this method is well known in gamma- 
gamma physics [7]. 

The main point is that the beams of virtual H^-boson are described not only by proba- 
bilities to find VF-boson with longitudinal and transversal polarization but by density matrix 
(as it should be in quantum mechanics). We consider the nondiagonal terms in the virtual 
W-boson density matrix, which turn out to be of the same order of magnitude as the diag- 
onal ones (corresponding to the longitudinal and transversal polarizations). If e + e~ (or qq) 
lab. frame and W + W~(ZZ) s.m.f. are collinear the contribution from these nondiagonal 
elements of density matrix to the total cross section is equal to zero. This is the reason 
why nondiagonal terms were not considered preveously. But e + e~ (or qq ) lab. frame and 
W + W~(ZZ) c.m.f. are not exactly collinear and the integration over azimuthal angle in the 
lab. frame does not cancel the contribution from non-diagonal density matrix terms. This 
means that two W boson beams can not be considered as an independent, i.e. factorization 
property is violated in this kinematical region. We have found that these interference terms, 
which having a negative sing and a large absolute value, explain the discrepancy between 
EWA and exact computer calculation. To restore factorization (and EWA method) we have 
to confine kinematical region to smaller scattering angles. We demonstrate that the EWA- 
method correctly applied yields appropriate lower bounds for the total cross-section for any 
energy of colliding particles and any (small or large) mass of Higgs-boson. The contribution 
of each element into the cross-section depends on the dynamics of the WW — > "X" subpro- 
cess (where "X" can be a Higgs boson(s), a,W,Z- boson pair or a lepton pair). For example 
for the WW — > H process (described as an example in this article) the contribution of lon- 
gitudinal and spin-flip interference terms are the most important ones and the transversely 
polarized term's contribution is negligible (contrary to [4] P.Johnson et al.). 

The paper is organized in the following way: the helicity density matrix for virtual W 
and Z bosons are defined in Section II. To illustrate our general consideration we use simple 
and well-known example - a single Higgs-boson production (see Section III). In Section IV 
we have discussed the exact numerical calculations of the contribution of each term in the 
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helicity decomposition into the total cross-section. In Section V we define the kinematical 
region for the application of the effective I¥-boson method. The results are discussed in the 
Conclusions. 

Some results presented here, have been already published in our recent preprint [8]. In 
[6] we have used this approach for two Higgs boson production via WVF-fusion. 



II. VK-boson density matrix 

To describe the polarization properties of the virtual W-boson we introduce the basis of 
virtual W-boson helicity states. It could be done in the center of mass system (c.m.s.) of 
W + and W~ 

q 1 = (iui,0,0,g) ; q 2 = (w 2 , 0, 0, -q) . 

Consider the standard set of orthonormal four- vectors orthogonal to the momentum q^. 

=(1)^^ = ^(0,-1,-^0), 



>(!)/ 



(0) 



def 



-It 



f (g,0,0,wi) = -iQi , 



(1) 



75(0,1,-^,0) 



and the analogous set for the second IT^-boson 

e? ) (±l)=ei 1) (Tl), 



f(0) 



-Q 2 2 



-q,0,0,w) d = iQ 2 , 



(2) 



e ]I' ■ fry. = • 



These four- vectors obviously represent the ±1 and helicity states of virtual W^-bosons in 
their c.m.s. and form a complete orthonormal basis for the subspace orthogonal to q^ and 
q 2il respectively: 

-l) a 5 a , b , a, b = ±1,0, (3) 



< ' : [<i)( ),(b) 



(4) 



Any vector or matrix orthogonal to q lfl could be expanded over this basis. Further we shall 
neglect the masses of light quarks and leptons. In this approximation 



and as a result 



a,b 







a = ±1, , 



P% = (-l)^e£\a)[4 1 \b)]*pV 
Here p a b is the density matrix in the helicity representation: 



(5) 

(6) 
(7) 



Pab — < 



P++ 


P+o 


p+- 






Pol 


(i) 
Poo 


(0 
Po- 


> , 


i = 1, 2 


(*) 
P-+ 


o {i) 

P-0 


(0 







(8) 
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(i) 

Detailed expressions for Lorentz invariant density matrix elements p a l are presented in the 
Appendix I and II. This density matrix is nondiagonal. It means that the system is strongly 
polarized and quantum mechanical interference effects are large. The diagonal components 
p++, P-L and Poo represent the fraction of longitudinal and transverse ly-bosons inside 
the fermion (electron) (unnormalyzed classical probabilities). The nondiagonal components 
p+o, Po_ and p+L correspond to the spin-flip transitions of VT-boson and depend on the 
azimuthal angle <pi in the W + W~ c.m.s. For the relative azimuthal angle between the 
scattering planes of the colliding particles Atp = (p\ — ip 2 , we have: 

COS A(y? = 1 , (9) 

pl± ■ pI± 

where pi± (p2±) is the perpendicular component of the electron momentum. [] 

Vy-bosons are produced by the left fermion component only, that is why the W^-boson 
density matrice elements are different from the photon density matrice ones [7], as opposed 
to photons, which produced by the left and right fermions with the same weight, and to Z- 
bosons ones (because Z-boson is produced by left and right fermion components, but with 
the different weights). These tiny effects are interesting to study at polarized e + e~-beams. 
We shall discuss them in the next publication. 

III. Exact calculation of Higgs-boson production 

The process of Higgs-boson production via PF-boson fusion in e + e~ (or qq) collisions is 
shown in (Fig. 1). Its cross-section can be written in the form: 

f-Ka w \ 2 1 q\ q\ 

d(7 ={— ) • — = „ „ —2 —2 z^rr ( 10 ) 



2 ' 4y / (p 1 p 2 ) 2 -p\pl (Qi-m 2 J 2 (ql - m 2 J 2 



4p2pS1 M; a M v ^Y5\ qi + q 2 - k)- 



d 3 k d 3 p[ d?p' 2 



(2n) 3 2E k (2rr) 3 2E[ (2n) 3 2E' 2 ' 

where a w = sin " g (a - is the fine structure constant and 9 W is a Weinberg angle); M^ a - 
denotes the amplitude for Higgs-boson production in the virtual W-hoson collision, in the 
Standard model: 

M m = gm w ■ g m . (11) 
Note, that the tensor structure - g^ a of this vertex is fixed by the unitarity and is thus the 

^ and Pal 



same in all models, except for some scaling coefficients [9]. The tensors pW and p^l are 



defined by the currents of the colliding particles (electrons): 



P$ = -^ 2 Tr[{p l +m e ) ltl {p' 1 +m e ) lv {l + 1 ,)}= (12) 
zq 1 



' Qi^iu\ (2pi - ?i) M (2pi - qi)„ m 2 e 

9^iv 2 — J 2 ' le LivapPlaQip + 9v^T 

v Qi J Qi Qi 



1 Note, that the iy-density matrix is defined only by the fermion- W- fermion vertex and 
have no dependence on the Higgs-boson models. 
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and represent spin density matrix for effective lU-bosons. 

If we substitute (11, 12) in formula (10) we shall obtain the known [10, 11, 13] expression 
for the matrix element squared: 



P^pS] M; a M uP = il(p lP2 )(p' lP ' 2 ) (13) 



where Pi,P2 (p>i,ph) are the four-momenta of the incident (final) electrons (quarks) (Fig.l), 
we have verified that the exact formula (10) with (13) gives numerically the same result 
[12] as the previous calculations [2, 10, 11]. (One integration with 5-function has been done 
analytically, see Appendix III.) Below, we shall discuss the contribution of the different 
W^-boson polarization states in the total cross-section. 

IV. Exact calculation of the input of different 
IV-boson polarization states into the total cross-section 

Using the helicity representation for pfy and pffi (6) - (8) one can rewrite the total 
cross-section (10) in the following way [6] 



1 (na w y q{ q$ r— r (1) (2) , 



i 4_ 

4(piP 2 ) V 2 J {q\ - m 2 w ) 2 {q\ - m 2 w f 

H- (tfjll + ) (Jrprp + 2 ( p W p S + h . C ) J LT+ (14) 

+ (pU p?> + h.C.) J TT ) ^d^dq\dqf^- , 
where ( = i] = A0 = X — 2 ~ is the relative azimuthal angle in e + e _ - c.m.s. 

a , (<?1-L • Q2±) / 1t ,x 

cos A<p = , (15) 

where gij^ffej.) ~ is the perpendicular momentum of IU + (W^) in e + e~ - c.m.s., \[X - is the 
flux of the virtual W^-bosons: 

X = (gig 2 ) 2 - q\ql , (16) 

a LL , and a TT represent the "cross-section" of the corresponding virtual lU-boson components 
[6]: 

n 2 a w 2 (qiq 2 ) 2 A/f2 \ /-^ 

GLL = ^/x mw HqT ' 4 n s - M *) > ( 17 ) 

^tt = • 45 (S - M|) (18) 

and J7" L t , Jtt ~ represent the interference terms [6] : 

JiT = ^,^|£.«(»-Mj), (19) 

J7rr — &TT , (20) 
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we have verified that exact calculation of the total cross-section (14) with the substitution 
(15-20) gives exactly the same result as formula (10) [12] (see Appendix IV). 

Now, being sure of the normalization, we can look at the contribution of each (17-20) 
term, multiplied by the corresponding density matrix element in (14), into the total cross- 
section. The numerical calculation of the contribution of the longitudinally and transversely 
polarized VT-bosons into the total cross-section (14) (a LL + a TT ) exact / 'cr^f is plotted in 
Fig.2 for M H = 100 GeV as a function of y/s. One can see that it is systematically larger 
than the total cross-section. This effect is easy to understand if we look at Fig.3, where the 
contribution of each term in (14) is shown, along with the total cross-section. The main 
contribution comes from the longitudinally polarized W^-bosons, while the contribution of 
the transversely polarized W-bosons is small. Interestingly, the spin-flip term Jlt, provides 
a large contribution with a negative sign which cancels a large part of the longitudinal 
contribution of a LL into the a to tai- 

The spin-flip terms p+o and p+_ in (14) depend on A0 - the relative azimuthal angle 
defined in the W + W~ c.m.s. (9), which is different from A0 (15) in the phase space volume 
(14) (defined in the e + e~ c.m.s.). 

Would coincide with A0, the contribution of interference terms into the total cross 
section was exactly zero. But these two azymuthal angles coincide only if two cm. systems 
moves along their Z axis. In general case A0 depends on A0 (see Appendix IV) and the 
contribution of quantum mechanical interference terms into cross section is non zero and 
rather large in the region of not very high energy and small Higgs-boson mass (Mh ~ mw) 
[12,13]. This subtle point was missed in the literature. 

From exact formula of Appendix III, taking 

U = ^9i9 2 (l - cos A(/?) + 1(0! - 6 2 ) 2 (21) 

we estimate, that e + e~ lab. system and WW-c.m.s. becomes acollinear when fermion 
scattering angles 9 li2 are greater then: 

4m ^ 

K > -jr^ ■ (22) 

s I 1 iL 



The nonvanishing contribution from interference terms J int in the total cross-section, com- 
pared to the longitudinal polarization contribution will be then: 

a LL M 2 H lns/M 2 H ' 1 ' 

i.e. it vanishes both at small scattering angles, and when M H ^> m w . 

V. Effective vector— boson method 

The standard version of the effective W-boson method [1, 2] could be derived from (14) 
if one takes into account only the terms proportional to a LL and a TT . Since (14) is already 
written in the factorized form, it is easy to calculate the effective FT-boson distribution 
functions. Of course they coincide with the standard ones. It is important now, to apply 
these functions only in the kinematical region, where the new additional terms in (14) (which 
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represents quantum mechanical interference effects) are not important. From Section 4, we 
have found that it can occur in the case of small scattering angle of electrons (22). It means 
that there must be a cut off in the integral over dqf, as q\ is connected with cos Of. 



gi 2 = --(l-cos^)-(l-C). 



(24) 



We have performed computer calculations for each helicity component in the total cross- 
section a tot and for a tot itself for different value of q 2 . The result is presented in Fig.4. 
We see that the transversal cross section ott is small compared to o tot for any q 2 , and 
its contribution can be neglected. In the region q 2 > Ml = 100 GeV 2 interference term 
-TLT{q 2 ) coincides with and the sum (<Jll + &TT){q 2 ) doesn't correspond to the o tot for 
this value of q 2 as it was always considered in the previous calculations. So the natural 
cut-off for q 2 is A 2 = s, where s is a characteristic energy of W + W~ cross-section. The same 
restriction on q 2 is needed to have W-boson beams to be independent from each other [6]. 
(For effective photons this cut-off was considered much earlier in [7].) 

For q 2 < s = Ml we have the following approximate formula for ull and <jtt'- 



appr 



a 



sin 2 9„ 



i+ 



M 



H 



In 



M 2 H 



-2 1- 



M 2 H 



(25) 



appr 
TT 



a 



m: 



sm 2 w J lQm 2 w M H 



2 1 + 



2M 2 H 



+ 



2s 2 



In 



M 2 H 



-3 1- 



M* H 




A 2 



(26) 



If we'll absolutely neglect (—tlx) i n this region q 2 < A 2 then we can calculate approximate 
total cross section. The ratios of ^TV^totaf an d (°ll + a TT) appr / &totai as functions of ^/s 



are plotted in Fig.5(a,b). We have taken cut-off exactly A 2 = Ml) and reproduce af^f 
with rather good accuracy. In principle one can also analytically take into account the small 
contribution of interference term — tlt (23). In this case we expect to obtain a little bit 
smaller cross section just because <r^ pr is obtained from a small part of the phase space, in 
the same way as we already know from 77-physics [7, 14]. 

Note, that the same formulas (25-26) are well applied in the case of high energy and 
heavy Higgs M H 3> m w with A 2 = Ml, and give the same result as in [2, 3, 10, 11] with a 
strong suppression of the transversely polarizeed W-heams, as 



O'LL 



M; 



II 



(27) 



Conclusions 

We have demonstrated that the approximate effective W^-boson method can be well 
applied in the kinematical region where the W + W~ c.m.s. is collinear with the c.m.s. of 
the colliding particles, i.e. in the case of small scattering angles for the electrons (or quarks) 
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and gives the lower bound of the cross-section. It means that there must be introduced a 
cut-off A 2 < Mfj when M H 3> m w ) or A 2 < (when M H ~ rn w ) in the integrals over the 
q 2 — W^-boson virtualities. In this region the contribution of the transversely polarized beams 
of W^-bosons are always smaller than the leading contribution of longitudinally polarized W- 
beams. 

It was shown, that when these two c.m.s. are not the collinear ones, the negative sign 
contribution of interference terms plays a crucial role and must be taken into account. That 
fact was missed in the previous works on this subject. 

It is important that the nondiagonal structure of W^-boson density matrix is defined only 
by the fermion-W^-fermion vertex and does not depend on the details of W^W^-interaction. 
But of course the contribution of each term of the density matrix in the particular subprocess 
will depend on the dynamics of the H^W-interaction and its particular kinematics. So in 
each particular case all the terms (not only the transversely and longitudinally polarized 
ones) of the density matrix must be taken into consideration. 
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Appendix I 



The virtual W^-boson density matrix elements: 



0- 



£1 

X 

jL 

2X 



£2 _ (9192) 



u 2 



e 2 ^ s 2 \ 
2 X [ 



f -77 



<i\il 



" ■ _|_ si _|_ (gig;) 



w s 2 
V2 X 



s 



H*-(1-*)(i-*-1) 



1/2 



Appendix II 



The virtual VT-boson density matrix in the case of small scattering angles of the colliding 
fermions [10]: 











2 










p'll 


= -2v^ie^^(l - 
C 2 


C) 1/2 , 




= -2v / 2^ 1 ^(l - 


C) 3/2 , 




= -2e 2 ^^(l-C) 





where ( = y=-. The analogous relation for p^ fe can be found with a substitution -^j — ■> -^2 
77; — > — </? 2 (and g| — > g 2 in Appendix I). 
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Appendix III 



For the exact numerical calculation of the total cross-section (10-13), we used the follow- 
ing final formula: 

-iK^ = ^L dcos9i L d cos 62 L h/s dC • 
i (i-^) , 

(l + ^-cos^l + ^-coe*)' V + <l-V)C 
where f]o = m "^^ J ^ U ] U = \ (1 — cos 9\ ■ cos 8 2 — sin 9\ sin 8 2 cos A(p) 

(0 < U < 1). 



Appendix IV 

For the exact numerical calculations of the helicity components in the total cross-section 
(14) (with (9), (15)-(20) ) we used as the final formula: 



2-nda ml,at f l 



j dcos#i J d cos 9 2 j M 2 H d( 



dAcp 2 6 s 2 
I Q1Q2 

2(M) -^(|p«||pS| + 1^11^1)} , 



cos , 



where i] , U are the same as in the Appendix III; 

q 2 2 = ~(l-Vo)(l-cose 2 ); q\ = - Ul - 0(1 - cosflx); ( qi q 2 ) 



COS A(/9 = 



2^ 



^ f Ml _ «i _ £\ _ ^ f Mi + 1± _ 



'/0 .si .s .9 s . S 2 .,3 



(c-f)(c-^) + c-f] 1/2 [(»-f)(»-^) + »4 



1/2 



cos 2 A</5 = 2 cos 2 A<p — 1 
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FIGURE CAPTIONS 



Fig. 1 Feynman graph for Higgs-boson production via W + W fusion in e + e 
(or qq) collisions. 

Fig. 2 The exact numerical calculation of the contribution of only diagonal 
terms of VU-density matrix into the total cross-section. It is seen that 
without interference terms it overstimates almost 3 times the total cross- 
section. 

Fig. 3 The exact numerical calculation of the contribution of each polarization 
of the VU-amplitude in the total cross- sect ion. It is very well seen the 
cancelation of by J7lt - the interference term. 

Fig. 4 The exact computer calculation of the dependence of each polarization 
of lU-amplitude on the PU-boson virtuality q 2 . One can see that at 
q 2 = A 2 = Mfj there is the total cancelation between longitudinal and 
interference terms. It means that the integration over dq 2 in the approx- 
imate calculations via effective TU-boson method has a cut-off A ~ M H - 

Fig. 5 The ratio of the approximate effective PU-boson method calculation of 
al P L pr (a) and (off + of ™)<w (b) with q 2 < A 2 = M 2 H over the exact 
computer calculation a exact . The ratio is always less than 1 above the 
threshold, it means that the effective VU-boson method gives the lower 
bound of the cross section. 
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